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Quasiparticle excitations adversely affect the performance of superconducting devices in a wide 
range of applications. They limit the sensitivity of photon detectors in astronomy [1, 2], the accuracy 
of current sources in metrology [3] , the cooling power of micro-refrigerators [4] , and could break the 
topological protection of Majorana qubits [5]. In superconducting circuits for quantum information 
processing, tunneling of quasiparticles across Josephson junctions constitutes a decoherence mech- 
anism [6 9]. As relaxation and pure dephasing times of transmon-typo charge qubits now reach 
100 /is, understanding whether quasiparticle tunneling may already bottleneck coherence is of high 
interest. We integrate recent advances in qubit readout [10] and feedback control [11] in circuit 
quantum electrodynamics [12] to perform the first real-time observation of quasiparticle tunneling 
in a transmon qubit. We demonstrate quasiparticle-tunneling contributions to qubit relaxation and 
pure dephasing in the millisecond range. Thus, quasiparticle tunneling will not limit coherence for 
at least one order of magnitude beyond the state of the art. 



The preservation of charge parity (even or odd number 

of electrons) has historically been of primary concern in 
superconducting quantum information processing (QIP) . 
In the first superconducting qubit, termed Cooper pair 
box (CPB) [13], maintaining the parity in a small island 
connected to a reservoir via Josephson junctions is essen- 
tial to qubit operation. The qubit states jO) and |1) con- 
sist of symmetric superpositions of charge states of equal 
parity, brought into resonance by a controlled charge 
bias rig and split by the Josephson tunneling energy Ej 
(< Ec, the island Cooper-pair charging energy). Quasi- 
particle (QP) tunneling across the junction changes the 
island parity, "poisoning" the box until parity switches 
back or rig is offset by ±e [14]. QP poisoning in CPBs 
has been extensively studied, with most experiments [IS- 
IS] finding parity switching times of 10 /iS — 1 ms, and 
some > 1 s [19-21]. While these times are long com- 
pared to qubit gate operations (~ 10 ns), the sensitivity 
of the CPB qubit transition frequency wqi to background 
charge fluctuations limits the dephasing time to < 1 /iS, 
severely restricting the use of traditional CPBs in QIP. 

Engineering the CPB into the transmon regime Ej ^ 
Ec [22, 23] exponentially suppresses the sensitivity of 
Woi to charge-parity and background charge fluctuations. 
However, recent theory [6, 8, 9] predicts that QP tun- 
neling remains a relevant source of qubit relaxation and 
pure dephasing, particularly as improved understanding 
of dielectric loss [24] and the Purcell effect [25] has al- 
lowed reaching the 100 /US scale in cQED. To guide further 
improvements, it is imperative to precisely pinpoint the 
timescale for QP tunneling and its contribution to qubit 
decoherence. To date, only upper and lower bounds have 
been placed [23, 26]. 

Here, we transform a state-of-the-art single-junction 



transmon qubit into a real-time charge-parity detector, 
and demonstrate that QP tunneling does not preclude 
reaching the millisecond timescale in transmon qubit co- 
herence. Our qubit is controlled and measured in a 
3D circuit quantum electrodynamics (cQED) architec- 
ture [24], i.e., the same environment that the qubit ex- 
periences in QIP applications, without need for any ad- 
ditional electrometer or other circuitry. At the heart of 
our detection scheme is a very small but detectable par- 
ity dependence of the qubit transition frequency (up to 
0.04% of the average woi/27r = 4.387 GHz), obtained by 
choosing Ej /Ec — 25. 

Standard Ramsey fringe experiments provide the first 
evidence of QP tunneling across the qubit junction, as 
shown in Fig. 1 for a refrigerator temperature T,. = 
20 mK. Instead of the usual decaying sinusoid, we ob- 
serve two. Repeated Ramsey experiments always reveal 
two frequencies, fluctuating symmetrically about the av- 
erage woi (Fig. Ic). The double frequency pattern results 
from the qubit frequency sensitivity to charge, with QP 
tunneling events shifting the energy levels by ±e in the 
rig axis. The fluctuation in the frequency difference A/ is 
due to background charge motion slow compared to QP 
tunneling. The observation of two frequencies in every 
experiment shows that QP tunneling is fast compared 
to the averaging time (^ 15 s), but slow compared to 
the maximum 1/2 A/ ~ 5 /iS [27]. From the similar am- 
plitude of the sinusoids, we infer that the two parities 
are equally likely. Clearly, these time-averaged mciasure- 
ments only loosely bound the timescale for QP tunneling, 
similarly to Refs. 23 and 26. 

In order to accurately pinpoint the timescale for QP 

tunneling, we have devised a scheme to monitor the 
charge parity in real time (Fig. 2a), taking advantage of 
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FIG. 1. Bistability and drift of the qubit transition frequency, a, Ramsey fringe experiment (dots) and best-fit sum 
of two decaying sinusoids (curve). Tlie reference oscillator is detuned 1 MHz from tiie average qubit transition frequency 
ojoi/Stt — 4.387 GHz. b, Sketch of the charge dispersion of the first two levels of the transmon qubit, showing 2e periodicity. 
QP tunneling across the junction shifts rig precisely by e, resulting in two transition frequencies /e and fo (not to scale), c, 
Repeated Ramsey experiments (15 s each) show a symmetric drift of /e and fo around iJoi/Svr, arising from background charge 
motion. The frequency difference 2 A/ = fe ~ fo ranges from to 1.76 MHz (see also the Supplementary Methods). 



recent developments in high-fidelity nondemolition read- 
out [10] and feedback control [11]. Starting from |0), the 
qubit is prepared in the superposition state (|0)-|-|l))/-\/2 
with a 7r/2 y pulse at Wqi- The qubit then acquires 
a phase ±7r/2 during a chosen idle time At = 1/4A/, 
where the + (— ) sign corresponds to even (odd) parity. A 
second Tr/2 x pulse completes the mapping of parity into 
a qubit basis state, even — >■ |0), odd — >■ A following 
projective qubit measurement [10] ideally matches the re- 
sult A/ — 1 (—1) to even (odd) parity. Feedback-based 
reset [11] reinitializes the qubit to |0) and allows repeat- 
ing this sequence every Aioxp = 6 /iS. We note that this 
scheme realizes the charge-parity meter in the envisioned 
top-transmon architecture [28], in which a transmon is 
used to manipulate and readout Majorana qubits. 

The time evolution of charge parity is encoded in the 
series of results M (Fig 2b) . The time series has zero av- 
erage, confirming that the two charge parities are equally 
probable. Both the QP dynamics and the detection 
infidelity determine the distribution of dwell times ti 
and t-i (Fig. 2d). The measured identical histograms 
match a numerical simulation of a symmetric random 
telegraph signal (RTS) with transition rate Frts, masked 
by uncorrelated detection errors occurring with proba- 
bility (1 — F)/2. These two noise processes contribute 
distinct signatures to the spectral density of M (Fig. 2c). 
The best fit of the form 



shows excellent agreement, giving l/Fj-ts = 0.79 ms and 
F = 0.92. 

While the above scheme detects a characteristic time 
for QP tunneling, additional experiments are needed to 
distinguish QP tunneling events that cause qubit transi- 
tions from those that do not. For this purpose, we model 
the system with four levels \kl) {k and / denote the qubit 
and parity state, respectively), and rates T'^f., connecting 
them, with k {k') and I {I') the initial (final) qubit and 
parity state, respectively (Fig. 3b). Based on the iden- 
tical distributions of dwell times, we consider symmetric 
rates F- . 

To extract the rates, we measure the autocorrelation 
function of charge parity, conditioned on specific initial 
and final qubit states (Fig. 3). Conditioning on a first 
charge-parity measurement Mi = +1 postselects the 
qubit in |0) and even parity. A second measurement 
M2 follows a waiting time r. By conditioning also on 
A/2 ~ +1, we ensure that the qubit both starts and ends 
in |0). A second parity measurement, ending with M^, 
completes the sequence. The average result, once cor- 
rected for detector infidelity (see Methods), is the parity 
autocorrelation Rqq{t) — (P(0)P(r))oo, with first (sec- 
ond) subscript indicating initial (final) qubit state. Ne- 
glecting qubit excitation, i.e., setting Fqi = roi+rg° — 0, 
Roo{t) simply decays as exp(— 2rQgi). The exact solu- 
tion shows that this remains a valid approximation when 
including the measured Fqi = 1/6 ms~^, since the proba- 
bility of multiple qubit transitions in r is negligible. Sim- 
ilarly, we measure the parity autocorrelation with qubit 



3 




FIG. 2. Real-time measurement of QP tunneling, a, Ramsey-type sequence converting the qubit into a cliarge-parity 
detector. Tlie sequence (see main text) is equivalent to a vr rotation conditioned on odd parity, b, Snapshots of a typical 
measurement trace, 48 ms long. Dots are repetitions of the experiment in a, at 6 ^s interval, c, Double-sided power spectral 
density of M , obtained by averaging the squared Fourier transform of 45 consecutive repetitions of b. The best fit of equation (1) 
gives l/Frts = 0.794 ± 0.005 ms. Repeated experiments have a standard deviation of 0.09 ms. d, Histograms of dwell times 
for M = ±1. The grey curve is a model of RTS with symmetric rate Frts and detection fidelity F, extracted from the fit in 
the main panel, e, Optical image of a qubit with identical geometry [10] to that used in this experiment. See Supplementary 
Methods for several control experiments testing the measurement protocol. 



initially and finally in |1), Ruij) « exp(-2rf?t). To 
do this, we use the same conditioning, but apply a tt 
pulse after Mi and before Af2. Exponential decay fits 
give = 0.92 ±0.04 ms and l/TH = 0.70 ±0.06 ms. 

To quantify the contribution of QP tunneling to the 
measured qubit relaxation time Ti = l/Tio = 0.14 ms, 
we apply the same method, but condition on initial state 
|1) and final state |0). The ratio of QP-induced to total 
relaxation rates a = Ffg/Fio (Fio — Ff^ + Pfp) can be 
extracted from Rio{t — )■ 0) = 1 — 2a. The best fit of the 
model Riq{t) to the data, with a as only free parameter, 
gives 1/Ffg = 3.3 ± 1.0 ms and 1/F^g = 0.14 ± 0.06 ms. 
This result clearly demonstrates that QP tunneling does 
not dominate qubit relaxation at Tr = 20 mK, contribut- 
ing only 5% of qubit relaxation events. 

To facilitate comparison to theory, we perform the 
above experiments at elevated (Fig. 4). We observe 
that Fj-ts, Poo, rf°, and Ffg have similar magnitude and 
jointly increase with in the range 20 — 170 mK. How- 
ever, Ti remains insensitive to Tr until 150 mK. The 
observed sign reversal in Rio{t — >■ 0) near this tempera- 
ture (Fig. 4b) indicates that QP tunneling becomes the 
dominant relaxation process. 

The effect of QP dynamics on the qubit degree of 
freedom in superconducting circuits has been extensively 
studied theoretically [6, 8, 9, 14]. For transmon qubits, 
the predicted QP-induced relaxation rate is [6, 8] 

TTo^^V^AuJoi, (2) 

TT 

where Xqp — 7iqp/2fc'oA is the QP density n^p nor- 



malized to the Cooper-pair density, with vq = 1.2 x 
10 fim~ -3/zeV"^ the single-spin density of states at the 
Fermi energy [17] and A the Al superconducting gap. 
This relation holds for any energy distribution of QPs. 
For Tj. > 150 mK, the data closely match equation (2) 
using the thermal equilibrium Xqp = ^eT^I^'' 
and A = 170 ^eV, the value estimated from the normal- 
state resistance of the junction (see Methods). The sup- 
pression of Ff g at lower is much weaker than expected 
from a thermal QP distribution. Using equation (2), we 
estimate riqp — 0.04 ±0.01 [imT^ at — 20 mK, match- 
ing the lowest value reported for Al in a Cooper-pair 
transistor for use in metrology [29]. Improved shielding 
against infrared radiation [30] could further decrease riqp 
at low Tr, and will be pursued in future work. 

QP tunneling events that do not induce qubit transi- 
tions contribute to pure qubit dephasing. Calculations 
based on Refs. 6 and 9 predict F^^ ~ Ffg, in good agree- 
ment with the data (Fig. 4c). It is presently not un- 
derstood whether these QP tunneling events completely 
destroy qubit superposition states (case A) or simply 
change the qubit precession frequency (case B). In either 
case, in the regime of strongly coupled RTS valid for our 
experiment (Fgg, F^'j' < A/ [27]) the QP-induced dephas- 
ing time is 2/(Fgg + Ff^j") ~ 0.8 ms. For case B, this time 
would further increase in the weak-coupling regime (at- 
tained at Ej/Ec > 60) due to motional averaging [27]. 

In conclusion, we have measured the characteristic 
times of QP tunneling across the single junction of a 3D 
transmon by converting the qubit into a real-time charge- 
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FIG. 3. Rates of QP tunneling with and without qubit 
transitions, a, Pulse sequence measuring the autocorrela- 
tion function of charge parity. Two parity measurements M\ 
and M3 (see Fig. 2) are separated by a waiting time r. Post- 
selection on Ml = 1 [10] prepares the state |0e). Similarly, a 
measurement M2 = 1 at the end of r ensures that the final 
qubit state is |0). Afa will coincide with M\ only if the parity 
is unchanged. Inserting tt rotations after M\ and/or before 
A/2 allows measuring the parity autocorrelation for different 
combinations of qubit states. A preliminary measurement Afo 
(not shown) initializes the qubit in |0) by postselection. b, 
Diagram of the four energy levels with the modeled transi- 
tion rates (not to scale), c, Charge-parity autocorrelation 
Rhk'ij) for qubit in state |0) (dots), |1) (squares), or having 
relaxed from |1) to |0) (diamonds) during r. The average of 
the conditioned Afs is corrected for detection infidelity (see 
Methods). Fitting the solution of the rate equations, con- 
ditioned on initial and final qubit state, gives the inverse 
rates: l/Fog = 0.92 ± 0.04 ms, l/Ffj" = 0.70 ± 0.06 ms, 
l/Ffg = 0.14 ± 0.06 ms, l/Ffg = 3.3 ± 1.0 ms. 
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FIG. 4. Temperature dependence of QP tunneling 
times. a,b, Charge-parity autocorrelation functions i?oo(''") 
and i?io(r) at 20,80 (a), 150 and 170 mK (b). i?io(r 0) 
progressively decreases, indicating an increasing contribu- 
tion of QP tunneling to qubit relaxation, c, QP tunnel- 
ing times for the ground-state (l/Foo, dots) and excited- 
state (l/rf°, squares) manifold extracted from Roo{t) and 
i?ii(r) (not shown). Fits to Sa/ (/) (Supplementary Methods, 
Fig. S5) give similar values for l/Fi-ts (diamonds). Dashed 
curve: theory for Fqq [6, 9] for thermally distributed QPs 
and A = 170 fieV. d, Relaxation times with (l/Ffo, upward 
triangles) and without (l/Ffg, downward triangles) QP tun- 
neling, obtained from Rio{t) and the overall relaxation time 
Ti = l/Fio (squares). Dashed curve: equation (2) for thermal 
equilibrium. Error bars are 1 s.d. 



parity detector. First, probing charge parity every 6 /iS 
with a Ramsey-like sequence reveals a symmetric RTS 
with 0.8 ms characteristic switching time. Second, mea- 
suring the charge-parity autocorrelation function, condi- 
tioned on specific initial and final qubit states, distin- 
guishes QP tunneling that induces qubit relaxation from 
that which does not. We have shown that QP tunneling 
is not the dominant relaxation mechanism, contributing 
just 5% of relaxation events in our state-of-the-art trans- 
mon. Reaching the millisecond horizon in coherence will 
facilitate the realization of fault-tolerant QIP with super- 
conducting circuits. 

METHODS 

Device parameters. The transmon has Josephson en- 
ergy Ej = 8.442 GHz and charging energy Ec — 



0.334 GHz. Using the Ambegaokar-Baratoff relation 
EjRn = A/8e^ and the measured room-temperature re- 
sistance i?n, 300K — 15.2 kfl of the single Josephson junc- 
tion, we estimate A = 170 ^eV. The qubit couples to the 
fundamental mode of the cavity uJr/2T: — 6.551 GHz (de- 
cay rate k/2-k = 720 kHz) with strength g/2'K = 66 MHz, 
inducing a dispersive shift 2x/27r = —1.0 MHz. The 
qubit relaxation time Ti may be limited by the Purcell 
effect [25]. A simple estimate including only the funda- 
mental cavity mode gives 240 jjLS. The dephasing time, 
r2 = 10 — 25 /iS, is limited by background charge fluctu- 
ations (see Supplementary Methods). 

Experimental setup. Projective readout with 99% fi- 
delity is achieved by homodyne detection with a 400 ns 
pulse at Wr- — Xi aided by a Josephson parametric am- 
plifier [10]. The qubit reset is implemented with a 
home-built feedback controller based on a complex pro- 
grammable logic device (Altera MAX V) that integrates 
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the last 200 ns of the readout signal and conditionally 
triggers a tt pulse (Gaussian, a = 8 ns, 32 ns long) 2 /xs 
after the end of the measurement [11] (see the Supple- 
mentary Methods for more details). 

Extraction of QP tunneling rates. To convert 
(M3(T))fefc' into Ri~k'{T), we correct for the overall de- 
tection errors, distributed among readout (< 1%) and 
reset (~1%) infidelities, suboptimal At (< 2%), and de- 
phasing during At (remaining 1 — 3%). For this correc- 
tion, we first fit an exponential decay to {M3(t))qq and 
(M3(r))ii. The average of the best-fit value at r = is 
used to renormalize the data in Figs. 3c and 4a,b. The fit- 
ted decay times are l/2rQQ and l/2r^^, respectively. To 
extract F^g and Ffg, we fit the solution of equation (2) 
to Rio{t), using Ffg -|- Ffg = Fiq. Fiq is obtained from 
the equilibration time Tcq after inverting the steady-state 
populations -P|o),ss) -P|i),ss with a tt pulse: 
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|0>,: 



(-P|0),ss + -P|l>,ss 



(3) 



The total excitation 1 — -F'|o>,ss is obtained by measure- 
ment and postselection [11]. Equation (3) remains a 
valid approximation even for the highest temperatures in 
Fig. 4, when the populations of higher excited states be- 
come relevant. In this case, the populations f]o),ss) -P|i>,ss 
are estimated from the total excitation, assuming that 
the populations are thermally distributed [11]. Error bars 
for r^g, F^Q are calculated from the standard deviation of 
repeated Ti measurements and the fit uncertainty in a. 
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EXPERIMENTAL SETUP 

The device and the experimental setup are similar to 
those described in Refs. 1 and 2. Here we detail the 
changes wc made since these earlier reports. To lower 
the transmon temperature, we replaced the Al cavity 
with a Cu cavity [3] , improved thermal anchoring to the 
mixing chamber plate, and added low-pass filters (K&L 
Microwave 6L250-8000/T18000-O/O) on the input and 
output ports of the cavity. As a result, the transmon 
temperature decreased from 127 to 55 mK, correspond- 
ing to a reduction of total steady-state excitation from 
~ 16 to 2%, respectively. Because these changes were 
made simultaneously, wc cannot pinpoint the individual 
contributions to the improved thermalization. 

To perform qubit reset faster, we replaced the ADwin 
processor with a home-built feedback controller based 
on a complex programmable logic device (CPLD). The 
CPLD allows a response time, from the end of signal 
integration to the 7r-pulse trigger, of 0.11 //s. The to- 
tal loop time, from the start of the measurement pulse 
to the end of the triggered tt pulse at the cavity input, 
is 0.98 /is, a substantial improvement over the previous 
2.62 jjLS. However, a delay is added to reach 2 jis (~ 10/k) 
between measurement and conditioned tt pulse, ensuring 
that the cavity is devoid of readout photons. 

VALIDATION OF THE CHARGE-PARITY 
DETECTOR 

We here perform several control experiments to val- 
idate the use of the qubit as a charge-parity detector. 
First, the parity-to-qubit state conversion is tested with 
suboptimal choices of the Ramsey interval Ai (Fig. SI). 
As expected from equation (1) in the main text, the white 
noise level in Sm increases at the expense of the signal 
contrast as At deviates from the optimal choice 1/4A/. 
Remarkably, the extracted rate Tpts is approximately con- 
stant down to F ~ 0.4 (Fig. Sic). This is consistent with 
the model of charge parity as a symmetric RTS, with time 
constant determined solely by QP tunneling. 

In a second test, we replaced the Ramsey-like sequence 
with a single pulse, with rotation angle 6. Time series of 
M for 61 = 0, TT and it/2 are shown in Fig. S2a. The 



very high occurrence (^ 99%) of 1 (—1) for = (tt) 
equals the efficiency of reset, following each measurement 
M. For 9 = it/2, the qubit is repeatedly prepared in an 
equal superposition of |0) and |1), and the measurement 
produces uncorrelated projection noise. The spectra of 
these control experiments are compared to the QP tun- 
neling measurement in Fig. S2b, clearly showing that the 
observed RTS is due to the signal acquired during At. 

As a final test of the charge-parity detector, we subject 
the qubit to an externally generated RTS, similarly to 
Ref. 4. Symmetric RTS sequences with switching rate 
Ftt are generated in Lab VIEW and sent to an ADwin 
controller. The ADwin samples the RTS at 9 /zs interval. 
When the signal is -1-1, the ADwin triggers an AWG520 
(also used for reset [2]), which then applies a it pulse on 
the qubit. As a result, the measured qubit state in M 
is conditioned on the RTS state, mimicking the parity- 
controlled IT pulse implemented in Fig. 2. In all cases, 
the fitted rates Ffit match the programmed F^ within 
3% (Fig. S3). 



SLOW FREQUENCY FLUCTUATIONS 

The slow drift of the two frequencies evident in Fig. Ic 
is attributed to rearrangement of background charges on 
the sapphire substrate, which capacitively couple to the 
qubit islands. In support of this hypothesis, T2* fol- 
lows an upward trend with increasing A/ (Fig. S4a), 
with maximum (20 — 25 /is) at the charge sweet spot 
(A/,„ax = 880 kHz). Histograms of A/ (Fig. S4b) are 
in good agreement with a sampled sinusoidal function, 
as expected for charge-modulated qubit frequency. The 
measured is consistent in order of magnitude with 
that estimated from the frequency power spectrum, pro- 
portional to 1//^ '' over 2 decades. 
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FIG. SI. Intentionally imperfect detector. Power spectral densities similar to those in Fig. 2, but with Ramsey time At 
varied from the optimum Atopt = 1/4A/. The detection fidelity F (d) decreases for a, At < Afopt and b, At > Atopt (b), 
but the extracted switching time l/Frts remains approximately the same, as expected (c). Curves in c are obtained from a 
numerical simulation using l/Frts = 0.78 ms. 
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FIG. S2. Control experiment replacing the Ramsey sequence in Fig. 2 by a ^-rotation of the qubit. a, 1 ms 

snapshots of measurements taken at 6 fis interval, over 48 ms. From top to bottom: 8 — 0, tt, and 7r/2. For 6 = 0(7r), 
M — 1 (—1), with ~ 1% error, vastly improving on the fast reset demonstrated in Ref. 2. For 6 = tt/2, M is randomized by 
projection noise, b, Power spectral densities of a. For 9 = n/2 the spectrum equals the white noise level (dashed line) at the 
sampling time of 6 fis. For 9 — ot tt, the spectrum is an oder of magnitude lower than the white noise level. The pickup 
at 2.1 kHz is associated with the pulse tube of the cryogen-free dilution refrigerator, as it disappears when the pulsed-tube 
compressor is momentarily turned off. The charge-parity spectrum from Fig. 2d is added for comparison. 
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FIG. S3. Power spectra of software-generated sym- 
metric RTS, as detected by the qubit. Programmed se- 
quences of symmetric RTS with switching rate r,r are loaded 
into an ADwin controller that triggers a vr pulse on the qubit 
when the level is high (similar to Ref. 4). The qubit is then 
measured and reset to |0) every 9 fis. This sequence emulates 
a RTS process similar to QP tunneling. The detected and 
programmed switching rates match to better than 3% over 
the range 50 to 1000 fj,s. 
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FIG. S4. Qubit dephasing and charge noise, a, Scatter plot of T2 as a function of A/ (see Fig. 1). The upward trend 
with A/ suggests that background charge noise is the main source of qubit dephasing. b, Histograms of A/, obtained from 
30, 000 Ramsey experiments acquired at ~ 4 s interval, with 5 kHz binning. The distribution peaks at the charge sweet spot 
(A/max = 880 kHz), where the qubit transition frequency is least sensitive to background charge fluctuations, c, Frequency 
power spectrum obtained from data in (b). We fit the equation A/ f" to the data for / > lO"'^ Hz, giving best-fit values 
A — 8.1 ■ 10^ and a — 1.7 (solid line). Dashed lines indicate the slopes for a = 1,2. 
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FIG. S5. Charge- parity spectra at variable Tr- Equa- 
tion (1) in tiie main text is fit to eacli spectrum. To facilitate 
comparison, data and fits liave tlie wliite noise offset sub- 
tracted and are renormalized. The best-fit rates are plotted 
in Fig. 4c. 



